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INTERFEROMETRIC  DENSITY  MEASUREMENTS  IN  THE  ARC  OF  A 
PULSED  PLASMA  THRUSTER* 

K.  I.  Thomassen  and  D.  Tong 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts 


Abstract 

Spatial  and  temporal  measurements  of  the 
electron  density  were  made  using  a  He-Ne  laser. 

The  multiple  pass  Mach-Zehnder  interferometer 
detects  densities  above  10*5  particles/cm^. 

Using  these  data  and  free  streaming  solutions  to 
the  Boltzmann  equation,  which  model  the  expansion 
of  the  exhaust,  we  determine  the  ion  temperature 
and  total  charge  in  the  exhaust.  The  percentage 
ionization  can  then  be  found,  and  estimates  made 
of  the  gas  dynamic  portion  of  the  thrust. 

I .  Introduction 

Pulsed  electric  thrusters  have  recently  found 
use  in  satellite  control  application*  and  are 
increasingly  being  incorporated  into  design  plans 
for  various  satellite  missions.  To  increase 
their  use  and  versatility  requires  improvements 
in  performance,  and  hence  we  initiated  a  program 
of  diagnostic  measurements  and  design  changes 
on  the  Teflon  fuel  thrusters  used  on  the  Lincoln 
Ibqperimental  Satellite,  LES-62>3.  Results  of  the 
diagnostic  work  to  date  may  be  found  in  the  liter - 
ature^ ,5. 

Earlier  measurements  determined  electron  temper¬ 
ature  and  densities  in  the  exhaust,  downstream 
from  the  plasma  arc  which  forms  over  the  Teflon 
surface.  To  determine  the  gas  dynamic  pressure 
of  the  plasma  on  the  fuel  surface  requires  a  density 
measurement  in  the  arc.  A  calculation  of  this 
surface  pressure  was  needed  to  determine  the  total 
force  on  the  exhaust.  Circuit  measurements  of  the 
arc  current,  giving  the  surface  magnetic  field, 
showed  that  the  plasma  acceleration  could  be 
explained  by  the  magnetic  pressure  on  the  arc. 
Nonetheless,  there  is  sufficient  uncertainty  that 
an  estimate  of  the  gas  dynamic  contribution  was 
required  to  complete  the  model  of  this  thruster. 

In  this  paper  we  report  an  experiment  to  meas¬ 
ure  the  arc  density  with  a  He-Ne  laser,  and  give  a 
theoretical  model  for  the  expansion  of  the  exhaust. 
From  the  theory  and  experimental  data  we  determine 
electron  density,  total  electron  charge  in  the 
exhaust,  and  ion  temperature .  The  ions  control 
the  plasma  expansion  so  a  density  measurement  with 
distance  determines  their  temperature.  We  show 
that  the  electron  pressure  on  the  plasma  contributed 
to  no  more  than  a  third  of  the  total  force  on  the 
exhaust. 

II.  The  Experiment 

A  multipass  Mach-Zehnder  interferometer  as 
described  by  Rizzo^is  used  in  this  measurement  to 
increase  the  total  path  length  through  the  system 
and  improve  the  sensitivity.  The  physical  arrange¬ 
ment  of  the  mirrors  and  beam  splitters  around  the 


vacuum  chamber  is  shown  schematically  in  Fig.  1(a). 
Multiple -pass  mirrors  in  both  the  reference  and 
plasma  arms  equalize  the  path  lengths  in  these  two 
arms.  The  beams  recombine  at  an  angle  0  as  shown 
in  Fig.  1(b)  to  produce  fringes  which  slide  past 
the  pinhole  collimator  as  the  density  is  varied. 

The  entire  interferometer,  including  the  source 
and  the  photodetector,  was  mounted  on  a  solid  aluminum 
I-beam  for  rigidity.  Optical  components  with  micro¬ 
meter  adjustable  mounts  were  supported  by  a  3/4” 
thick  aluminum  yoke  bolted  to  the  I-beam.  Micrometer 
and  screw  adjustments  were  found  to  be  indispensable 
in  the  alignment  of  the  interferometer. 

Since  the  beam  passes  through  the  exhaust  plasma 
at  several  positions  along  the  arc,  we  obtain  an 
average  of  the  arc  density  (in  the  direction  of 
the  current) .  This  presents  no  problem  because 
we  are  interested  primarily  in  the  density  perpendic¬ 
ular  to  the  Teflon  face.  A  diverging  lens  is  used 
to  expand  the  fringe  pattern  before  it  reaches  the 
pinhole,  and  a  polaroid  lens  attenuates  the  fringe 
pattern  below  the  saturation  level  of  the  phototube. 

The  source  used  is  a  3mW,  c-w  He-Ne  laser 
operating  at  6328A.  The  beam  splitters  and  mirrors 
are  high  quality  1/20  X  optical  flats,  dielectric 
coated  to  operate  at  the  laser  wavelength.  The 
vacuum  windows  are  flat  to  1/10  X  and  have  anti- 
reflective  coatings  to  minimize  beam  attenuation. 

The  phototube,  Amperex  type  XP-1002,  has  a  response 
linear  to  within  one  percent.  Both  the  phototube 
and  the  oscilloscope  have  response  times  in  the 
nanosecond  range. 

In  spite  of  the  rigid  optical  specifications, 
the  fringe  pattern  degrades  to  an  unacceptable  level 
when  the  beam  is  reflected  for  more  than  7  passes 
through  the  plasma.  This  degradation  is  due  mainly 
to  the  vacuum  windows  and  not  to  the  coherence  length 
of  the  laser.  It  was  not  possible  to  include  the 
interferometer  inside  a  larger  vacuum  chamber  because 
the  exhaust  would  contaminate  the  optical  surfaces. 
Because  of  this  contamination  the  vacuum  windows 
had  to  be  cleaned  periodically. 

Vibration  of  the  interferometer  was  a  severe 
problem  when  the  vacuum  chamber  was  mounted  directly 
on  the  I-beam.  Vibrations  from  the  roughing  pump, 
transmitted  through  the  rubber  hose  connection  to 
the  vacuum  chamber  and  hence  to  the  interferometer, 
were  sufficient  to  destroy  any  observable  fringe 
pattern,  making  alignment  virtually  impossible. 
Eventually  the  interferometer  was  isolated  from 
the  vacuum  system,  reducing  vibration  to  an  acceptable 
level.  The  remainder  was  found  to  be  useful  in 
measuring  the  amplitude  of  the  fringes.  This  calibra¬ 
tion  was  needed  because  the  plasma  introduced  less 
than  2tt  radians  of  phase  shift.  No  errors  are 
introduced  in  the  measurement  by  the  vibration  because 
its  period  is  on  the  order  of  milliseconds. 


*"~kxperiments  were  performed  at  M.I.T.  Lincoln  Laboratory  and  sponsored  by  the  Department  of  the  Air  Force. 
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To  determine  the  electron  density  from  the 
detector  signal  we  combine  the  electric  fields 
of  the  two  waves  shown  in  Fig.  1(b), 
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where  k  *  2t r/X  is  the  wavenumber  of  the  laser 

source  and  a  and  3  are  the  z  and  y  compon¬ 
ents  of  k.  The  phase  difference  <f>  between  the 
two  paths  includes  the  plasma  contribution  which 
changes  by  A<f>  when  plasma  is  present: 


A<t>  =  kd(l  -  y  1  -  co  2/w2) 
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with  ojp  the  radian  plasma  frequency.  This  phase 
shift  is  proportional  to  density  in  our  experiment. 

The  width  d  must  be  replaced  by  an  effective 
width  if  the  density  varies  with  distance  along 
the  beam.  This  width  is  found  by  integrating  the 
density  profile  and  equating  it  to  the  peak  density 
times  an  effective  width.  The  fields  in  (1),  with 
their  associated  magnetic  fields,  give  an  average 
Poynting  flux  whose  z -component  is 

<pz>  =  hW +  e22  +  2EiE2(cos  ^ +  *1  (3) 

with  n  =  %/p0/e  =  377  ft.  The  fringes  along  the 

y  axis  and  their  shift  with  plasma  density  can 
be  determined  from  (2)  and  (3) ,  since  the  detected 
voltage  is  proportional  to  (3) . 

The  position  of  the  thruster  in  the  vacuum 
chamber  can  be  adjusted  by  an  external  connecting 
rod.  With  the  thruster  at  various  positions,  oscil¬ 
lographs  are  taken  of  the  phototube  output  by  trig¬ 
gering  the  thruster  and  oscilloscope  simultaneously. 
These  oscillographs  give  the  time  variation  of 
phase  shift  at  various  positions  downstream.  In 
Fig.  2  the  upper  trace  shows  intensity  variations 
due  to  vibration  while  the  lower  trace  shows  intensity 
variations  due  to  the  presence  of  plasma  (note 
the  different  time  scales) . 

The  effective  plasma  width  de  used  at  various 
downstream  positions  is  given  in  Table  I  from  the 
theoretical  results  of  the  next  section.  At  distances 
far  from  the  Teflon  face  the  arc  looks  like  a  point 
source  so  the  free  streaming  solution  is  a  simple 
gauss ian  in  the  coordinates.  (See  Eq.  7.)  Then, 
de  increases  linearly  with  time  and  is  given  by 


Table  I  Plasma  effective  width  used  in  the 
calculation  of  maximum  density 


Distance  from 
Teflon  Face 

(mm)  

de  (mm) 

1.0 

3.0 

3.5 

3.0 

6.0 

3.6 

8.5 

5.1 

11.0 

6.6 

13.5 

8.1 

16.0 

9.6 

21.0 

11.9 

26.0 

14.7 

By  the  methods  described  above  we  measured 
the  density  as  a  function  of  time  at  various  distances 
along  the  thrust  axis.  The  curves  given  in  Figs. 

3  and  4  are  averages  over  data  taken  from  several 
oscillographs.  In  each  case  the  beam  made  5 
passes  through  the  plasma.  An  average  is  required 
because  the  amplitudes  vary  from  discharge  to  dis¬ 
charge.  The  general  shape  of  the  curves,  however, 
remains  about  the  same. 

Maximum  electron  density  at  the  Teflon  face 
occurs  approximately  0.6  microseconds  after  initi¬ 
ation  of  the  discharge  when  the  current  is  a  maximum. 
Assuming  de  =  3mm,  the  peak  density  is  approximately 
9.3  x  lolfycc.  Peak  densities  downstream  are  less. 
Figure  5  shows  the  peak  electron  density  (in  time) 
at  any  downstream  point  as  a  function  of  distance 
downstream.  The  peak  density  decreases  slowly 
at  first,  but  eventually  varies  inversely  as  the 
cube  of  the  distance.  The  lowest  density  (at 
20  cm)  was  measured  with  the  microwave  interferom¬ 
eter  3.  Unfortunately  densities  between  26  mm  and 
20  cm  had  to  be  extrapolated  since  neither  inter¬ 
ferometer  operates  in  that  range. 

All  of  the  oscillographs  taken  with  the  beam 
adjacent  to  the  Teflon  face  show  a  reversal  of  phase 
shift  indicating  that  the  dielectric  constant  changes 
from  a  value  less  than  zQ  to  one  greater  than  e0. 
Such  a  reversal  is  not  seen  in  any  of  the  oscillo¬ 
graphs  taken  further  downstream,  so  neutral  atoms 
are  assumed  responsible.  If  this  is  the  case  an 
estimate  of  the  neutral  density  near  the  Teflon  face 
can  be  made. 

The  dielectric  constant  of  a  gas  of  neutral  par¬ 
ticles  is  given  by 


=  e  (1  +  — ) 

oy  Z  J 

o 


d  =  /iF  Vrp  t 

e  T 

where  Vj  is  the  ion  thermal  velocity.  At  distances 
on  the  order  of  the  arc  width,  d  is  approximated 
by  the  arc  width  of  3  mm,  as  determined  from 
high  speed  end-on  photos  of  the  plasma  luminosity. 


where  N  is  the  particle  density  and  a  is  the 
polarizability,  defined  as  the  ratio  of  the  induced 
dipole  moment  of  an  atom  to  the  applied  field.  In 
general,  a  includes  several  effects:  electronic 
polarizability  due  to  separation  of  the  electron 
cloud  from  the  neucleus ;  ionic  polarizability  due 
to  separation  of  atoms  in  a  molecule  with  ionic  bond¬ 
ing;  and  orientational  polarizability  due  to  the 
tendency  of  molecules  with  permanent  dipole  moment 
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to  orient  themselves  along  the  applied  field  direc¬ 
tion.  In  the  ablated  teflon  gas,  we  expect  mostly 
single  atoms  to  exist  because  of  the  low  molecular 
bonding  energy.  Hence  only  the  electronic  polariza¬ 
bility  of  carbon  and  fluorine  atoms  contributes 
to  the  dielectric  constant. 

A  very  simple  model  represents  the  electrons 
of  an  atom  by  a  sphere  of  uniform  negative  charge 
density  about  the  positive  neucleus?.  From  this 
model  one  finds  that  a  is  proportional  to  the 
atomic  volume.  Using  this  result,  the  polarizability 
for  the  actual  gas  is  predicted  to  be  3  x  lCT^l 
Farads -m^,  a  value  between  the  experimental  values 
for  helium  and  neon. 

The  magnitude  of  the  reverse  phase  shift  is 
0.2  rad.,  corresponding  to  e  =  eQ(l  +  2.6  x  10"^). 
Using  the  assumed  value  of  a.  we  calculate  a 
neutral  density  of  7.6  x  10^°/cc. 

III.  Expansion  of  the  Exhaust 


just  given  with  tne  source  function,  determined  by 
measurement.  These  solutions  can  then  be  fitted 
to  density  measurements  downstream  to  determine 
the  parameters  of  interest,  which  we  now  proceed 
to  do. 


For  a  given  group  of  particles,  the  equilibrium 
distribution  to  which  they  relax  via  collisions 
is  the  Maxwellian.  Further,  while  the  arc  is 
forming  the  particles  are  accelerating,  but  when 
the  arc  formation  is  complete  we  assume  the  par¬ 
ticles  have  acquired  a  constant  drift  velocity 
super imposed  on  the  Maxwellian  distribution. 

The  distribution  function  is  then  of  the  form 


f  «  exp 


1 


+ 


(8) 


where  the  ion  thermal  velocity  is  vj,  and  vQ 
is  the  drift  velocity  along  the  thrust  axis  (x 
axis) . 


The  plasma  formed  in  the  discharge  arc  expands , 
free  from  collisions,  into  the  vacuum  surrounding 
the  thruster.  This  expansion  is  governed  by  the 
expansion  of  the  ions,  hence  a  free  streaming  solution 
of  the  Boltzmann  equation  for  the  ion  distribution 
determines  the  parameters  describing  the  exhaust. 

In  particular,  comparing  our  measurements  with 
the  results  of  the  theory  gives  both  the  average 
ion  exhaust  velocity  and  ion  temperature. 

Before  giving  the  free  streaming  solutions 
we  show  that  a  collision  model  gives  answers  which 
are  in  contradiction  with  the  interferometer  data. 

The  force  free  Boltzmann  equation  with  a  simple 
collision  model  is 

II  +  V  •  Vf  =  -  v(f  -  fG)  (4) 

From  the  first  two  moments  of  this  equation  one 
obtains  the  usual  equations  of  continuity  of  mass 
and  motion,  which  combine  to  give 


We  are  interested  in  the  region  x  >  0  where, 
initially,  f  is  identically  zero.  Plasma  is 
being  generated  in  the  arc  which  is  placed  at 
x  <  0.  The  boundary  x  =  0  can  be  chosen  arbitrar¬ 
ily  as  some  point  where  density  data  r^ft)  is 
available.  With  the  assumption  made  above  on 
the  velocity  distribution,  the  boundary  condition 
on  f  is 


(9) 


where  g(y)  and  h(z)  are  as  yet  unspecified 
functions  and  A  is  a  normalizing  constant. 

This  expression  holds  only  for  vx  >  0  and  is 
identically  zero  for  v^  <  0.  Physically,  all 
particles  enter  the  region  x  >  0  through  the 
x  =  0  plane  with  positive  velocities.  This  restric¬ 
tion  becomes  important  when  f  is  later  integrated 
over  velocity  to  obtain  the  particle  density  n. 


+  DV2n  =  0  (5) 
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where  n  is  the  ion  density,  D  =  v~,  t  is  a  diffu¬ 
sion  coefficient  describing  the  collision  process 
(relaxation  time  x) .  Solutions  to  (5)  assuming 
an  initial  point  source  of  particles  is 


n  <*  t'3//2  exp 


r  x2  ♦  y2  ♦  z2l 

r  <7/ 1  x  J 


(6) 


Solutions  to  Eq.  4  in  the  absence  of  collisions 
are  obtained  by  the  method  of  characteristics.  For 
all  x  >  0,  f  is  given  by 


(10) 

as  can  be  verified  by  direct  substitution.  It  follows 
then  that  the  solution  we  seek  is 


The  free  streaming  solutions  to  (4)  (with 
v  =  0)  are 


f(r.v,t) 


An  (t 
ov 


V  x 

h(z  -  y5-) 
x 


_.-3  ("•  x2  +  y2  +  z2l 

n-t  “pr  m 

For  times  t  >  4i  the  free  streaming  solution 
dominates.  During  arc  formation  of  course,  collis¬ 
ions  dominate  and  establish  a  quasi -equilibrium. 
Consequently,  to  describe  the  exhaust  we  must  use 
the  free  streaming  solutions.  In  addition,  since 
the  plasma  is  created  at  the  surface  during  a  finite 
time  period  we  must  convolve  the  impulse  solution 


which  satisfies  the  boundary  condition  at  x  =  0. 
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Theoretically,  we  could  integrate  £  over 
velocity  to  obtain  the  density  as  a  function  of 
space  and  time,  but  since  the  exact  forms  of  g 
and  h  are  unknown  integration  is  impossible. 

A  different  approach  is  taken  here  whereby  the 
half  space  x  >  0  is  divided  into  three  sections; 
the  near,  mid,  and  far  regions.  In  the  near  region, 
adjacent  to  the  discharge  arc,  the  source  is  planar. 

In  the  mid  region  the  arc  appears  to  be  a  line  source, 
while  the  far  region  is  where  the  arc  is  a  point 
source.  Boundary  conditions  in  these  regions  are: 


noAi  exP[v| 

VS;  exP  PH  5(z) 

n0A3  exp  nn  6(y)  <$(z) 


(near  region) 
(mid  region) 

(far  region) 


where 


exp  [v]  =  exp 


-71  h  -  Vo> 


2  2  2 
+  v  +  v 
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In  each  region  f  is  still  given  by  Eq,  10. 


Since  nQ(t)  is  an  experimentally  determined 
function,  n(r,t)  can  be  obtained  only  by  numerical 
integration  of  f.  Alternatively,  it  is  equally 
valid  first  to  obtain  the  solution  to  a  time  impulsive 
source,  and  then  to  convolve  it  with  n  (t)  to 
obtain  n(r,t).  One  can  obtain  considerable  informa¬ 
tion  from  the  impulse  response  without  carrying  out 
the  convolution.  Therefore,  Table  II  lists  the 
distribution  function  f^  and  density  n^  for 
a  time  impulse  source  in  each  of  the  three  regions. 

In  each  region,  the  actual  density  n(r,t)  is  given 
by 


n(r,t) 


/ 


t 

n0CO  -  t)  dx 


(12) 


where  n.  is  the  time  impulsive  solution  for  that 
particular  region. 

For  the  parameters  of  this  plasma,  the  duration 
of  the  time  impulsive  response  for  the  first  25 
mm  in  the  exhaust  is  short  compared  to  the  duration 
of  the  source  (see  Fig.  6).  Consequently,  the  con¬ 
volution  will  give  a  maximum  density  when  the  impulse 
response  occurs  at  the  time  when  the  source  function 
is  a  maximum.  This  maximum  is  simply  the  area  under 
the  curve  of  the  impulse  response  at  the  distance 
in  question,  times  the  peak  source  density.  To 
evaluate  the  integrals  of  n^(t)  in  Table  II  is 
difficult,  hence  we  obtain  a  lower  bound  on  this 
area  by  computing  the  peak  value  in  time  of  the 
impulse  response.  This  technique  predicts  x*l, 
x~2,  and  x_3  dependences  of  the  peak  density  through 
the  three  regions.  For  the  far  region  this  technique 
is  exact  because  the  actual  density  approaches  the 
impulse  solution. 

The  mid  and  far  region  solutions  in  Table  II 
both  have  the  same  gauss ian  dependence  in  the  trans¬ 
verse  dimension.  The  interferometer  measurements 
give  no  information  on  the  transverse  variation 
of  density,  but  from  previous  Faraday  cup  current 
measurements  this  data  is  available.  The  maximum 


of  the  collected  current  (in  time)  as  a  function  of 
distance  across  the  beam  is  a  gaussian  with  a 
1/e  width  of  8.4  cm  at  a  distance  of  15.2 
cm  from  the  arc.  To  test  our  theoretical  results 
we  will  derive  an  expression  for  that  maximum 
current . 

The  ion  current  in  the  x  direction,  Jx, 
is  related  to  the  distribution  function  f  by  the 
equation 


■  ‘f/f 


vxf  dv 


Carrying  out  the  integral  using  the  far  region  dis¬ 
tribution  function  yields 


eA 


The  time  when  the  maximum  current  occurs  at  a  fixed 
position  x  is  given  by 


t  =  51 

mvJ  v 


where ,  for  M 


vo/vT 


MC 


M2 

r* 


fW$ 


(13) 


Far  from  the  source  but  near  the  thrust  axis  (z  «  x) 
the  peak  collected  ion  current  perpendicular  to 
the  arc  and  to  the  exhaust,  J^O^y  =  0,z),  is 
given  by 


Jmx(x*y 


0,z)  <*  x  ^  exp 


[-  (Z/C  x)2]  (14) 


This  gaussian  dependence  with  z  and  the  x 
dependence  at  z  =  0  was  verified  by  the  Faraday 
cup  current  measurement. 

An  estimate  of  total  charge  in  the  exhaust 
can  be  made  by  integrating  the  ion  flux  n(x,t)vx 
crossing  an  arbitrary  x -plane  to  infinite  time. 
Having  measured  n(t)  at  some  point  x,  we  can 
express  the  particle  distribution  function  at  that 
plane  using  Eq.  9, 

£  =  AT(t)  g(y)  h(z)  exp[v) 


where  T  is  proportional  to  the  measured  density 
variation  and  we  normalize  T,  g,  and  h  to  have 
a  peak  value  of  unity.  The  constant  A  is  related 
to  the  maximum  measured  density  n^.  Since 


/Aj  /eo 

dvx J  dvy J  dvz  fo 


T,g,h 


4 


we  have 


A  = 


rat 

7i  /^x)  /Jo  "  L  " 

dvxy  dvy  y  dvz  exp  [v] 


(15) 


-oo  -oo 


The  total  current  crossing  the  x-plane  is  given  by 

Jx  =  eAT(t)  g(y)  h(z)  j*  dvx  dvy J*  dvz  vxexp[v] 

(16) 


-OO  -oo 


Integrating  Jx  over  t,  y  and  z,  we  get  the 
total  charge  Q  passing  the  x-plane. 


//"bo  /3o 

dt  J  dyl  dz  J, 

-uu  -OO  -OO 


(17) 


Combining  Eqs.  15,  16,jand  17  and  carrying  out 
the  integration  over  v, 


Q =  e  Vvc 


,-M2 


1  + 


2H 


4  +  «P(M). 


//&  f oo 

dt  f  dy  /  dz  Tgh 

-oo  -oo  -oo 


where  M  =  vq/v^  and  4>(M)  =  exp(-s2)ds. 


IV.  Comparison  of  Theory  and  Data 

In  this  section  we  compare  our  calculations 
with  the  data  to  determine  the  applicability 
of  the  model  and  to  give  values  for  some  of  the 
exhaust  parameters.  The  predicted  variation  of 
maximum  density  with  distance  is  verified  and 
the  extent  of  the  three  solution  regions  is  found. 
Values  of  vp,  v0,  and  Q  are  then  calculated. 
Finally,  with  the  estimated  values  of  v0  and 
vr,  the  convolution  integral  given  by  Eq.  12  is 
carried  out  numerically. 

The  experimentally  determined  variation  of 
maximum  density  (in  time)  in  the  thrust  direction 
is  shown  in  Fig.  5,  where  the  curve  decreases  slowly 
at  first  but  ultimately  varies  as  the  inverse  cube 
of  the  distance.  From  the  plot,  using  the  expected 
x“l,  x"2,  and  x"3  dependences,  we  deduce  that 
the  near  region  extends  from  1  -  5  ran,  the  mid 
region  from  10  -  20  ran,  and  the  far  region  from 
25  ran  on.  By  comparison,  the  arc  dimensions  are 
about  25  ran  in  length  (y  direction)  and  3  ran 
in  width  (z  direction) . 

An  estimate  of  M  *  v0/vj  can  be  made  by 
examining  the  variation  of # peak  ion  current  across 
the  exhaust.  Comparing  the  Faraday  cup  data  with 
Eq.  14  gives  C  =  0.28.  Using  Eq.  13  with  z  «  x, 
we  find  M  ^  3. 


This  estimate  of  M  assumes  a  single  drift 
velocity  for  the  exhaust.  However,  there  are 
multiply  ionized  particles  in  the  plasma  with 
different  drift  velocities.  The  more  highly 
ionized  particles  drift  at  higher  velocities, 
and  since  ion  current  is  proportional  to  both  the 
particle  charge  and  velocity,  the  Faraday  cup  data 
is  weighed  in  favor  of  the  multiply  ionized  particles. 
This  weighing  is  evidenced  by  the  fact  that  the 
arrival  time  of  peak  current  corresponds  to  a  speed 
of  40  km/sec,  whereas  the  average  speed  of  all 
the  ionized  particles  is  approximately  30  km/sec. 

Hence  the  actual  value  of  M  is  probably  lower 
by  an  amount  dependent  on  the  percentage  of  multiply 
ionized  particles. 

Having  the  above  value  of  M,  we  can  derive 
vQ  and  vT  from  the  expression  in  Table  II  for 
tjjj,  the  arrival  time  of  peak  density.  Note  that 
in  the  far  region,  for  »1,  %  x/v0.  The 

maximum  density  at  20  cm,  as  measured  by  the 
microwave  interferometer,  occurs  at  6.5  microseconds. 
The  peak  arc  density  occurs  at  0.5  microseconds. 

Hence  an  elapsed  time  of  6.0  microseconds  is 
chosen  for  tm*  Using  the  far  region  expression 
for  with  M  =  3,  we  find 

vT  =  9.3  km/sec 

(18) 

vq  =  27.9  km/sec 

This  value  for  v0  is  in  good  agreement  with  experi¬ 
mental  evidence  so  far.  The  thermal  velocity  corres¬ 
ponds  to  an  ion  tenperature  of  7.6  eV  for  ions 
having  a  mass  of  16.7  atomic  units,  an  average 
for  2  carbon  and  4  fluorine  atoms.  The  estimate 
for  vj  is  low  since  the  estimate  of  M  is  expected 
to  be  high. 

With  these  values  of  v0  and  Vj,  the  convolu¬ 
tion  given  by  Eq.  12  can  be  carried  out  numerically 
to  obtain  predicted  density  as  a  function  of  time 
for  comparison  with  data.  Since  we  want  to  avoid 
the  region  where  the  arc  is  formed,  we  take  the 
left  boundary  plane  of  the  solution  region  to  be 
the  plane  at  6  mm  away  from  the  Teflon  face. 

The  solution  downstream  is  arrived  at  by  numerically 
convolving  the  density  data  at  the  boundary  with 
the  appropriate  time -impulsive  solution  as  given 
in  Table  II.  The  normalizing  constant  Ap(i  =  1, 

2,  or  3  corresponding  to  near,  mid,  or  far  region) 
is  found  by  integrating  f  over  velocity  and  matching 
the  result  to  nQ.  In  thecases  of  the  mid  and 
far  regions,  it  is  necessary  to  multiply  the  constant 
by  the  effective  plasma  cross  section  because  of 
the  spatial  impulse  approximation.  Hence, 


A1 

=  V3*3/2 

(near  region) 

*2 

=  deAl 

(mid  region) 

A3 

“  dedgAl 

(far  region) 

where  d^  is  the  effective  plasma  width  in  the 
z  direction  and  dg  is  the  effective  plasma  width 
in  the  y  direct  ion*  (parallel  to  arc) .  Only  the 
near  and  mid  region  solutions  apply  in  the  following 
calculations  since  data  is  not  available  more  than 
25  ran  from  the  Teflon  face.  Figure  6  shows  the 
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measured  density  at  6  mm  and  the  impulse  responses 
at  various  distances  downstream.  Calculated  and 
measured  density  curves  are  shown  in  Fig.  7.  The 
convolved  curve  at  8.5  mm  (Fig.  7a)  uses  the  near 
region  solution  while  the  convolved  curves  at  13.5 
mm  (Fig.  7b)  and  at  16  mm  (Fig.  7c)  use  the  mid 
region  solution.  In  each  case  the  curves  are  similar 
in  shape  but  the  theoretical  curves  peak  higher 
than  the  measured  ones.  This  discrepancy  may  reflect 
too  low  an  estimate  of  Vp,  the  neglect  of  recombin¬ 
ation,  or  that  the  source  was  approximated  by  an 
impulse  in  space.  Carrying  out  the  convolution 
with  vj  arbitrarily  increased  by  50%  did  not 
result  in  significant  change  in  this  difference. 

The  possibility  of  significant  recombination  is 
inferred  from  Figs.  7  where  the  ratio  of  measured 
peak  density  to  calculated  peak  density  decreases 
with  distance  from  the  Teflon  face. 

The  total  charge  Q  is  given  theoretically 
by  Eq.  18.  Since  we  need  to  know  only  the  areas 
under  the  curves  T,  g,  and  h  to  compute  Q,  we 
can  approximate  them  by  square  pulses  of  unit  magni¬ 
tude  and  some  estimated  effective  widths.  The  effec¬ 
tive  width  of  the  pulse  approximating  h(z) ,  for 
example,  is  the  value  of  de  given  earlier.  Since 
M  %  3,  the  quantity  in  brackets  in  Eq.  18  is 
approximately  1.  Hence  Eq.  18  reduces  to 


Q  *  eVo¥gde  (19) 

where  the  d's  represent  the  effective  widths. 
Estimated  values  of  Q  and  the  associated  effective 
widths  are  listed  in  Table  III.  For  comparison, 
the  total  charge  derived  by  integrating  over  the 
Faraday  cup  current  is  also  given.  The  steady 
decrease  in  estimated  values  of  Q  gives  further 
evidence  that  significant  recombination  occurs. 

V.  Summary 

An  extensive  study  has  been  done  at  M.I.T. 
Lincoln  Laboratories  on  the  solid  Teflon  plasma 
thruster,  including  the  development  of  a  circuit 
nodel  and  its  analysis,  and  several  plasma  diagnostic 
experiments.  However,  density  measurements,  using 
Langmuir  probes  and  a  microwave  interferometer, 
failed  to  reveal  the  particle  density  in  the  discharge 
arc,  without  which  it  is  impossible  to  estimate 
the  relative  importance  of  gas  dynamic  and  magnetic 
forces.  So,  the  laser  interferometer  described 
in  this  paper  was  constructed  and  used  to  make 
the  desired  electron  density  measurements.  The 
peak  electron  density  in  the  arc  was  found  to  be 
9.3  x  lOiG/cc. 

The  gas  dynamic  pressure  can  be  calculated 
assuming  the  ion  density  is  approximately  equal 
to  the  measured  electron  density.  We  determined 
that  the  ion  temperature  is  7.6  eV,  and  measurements 
of  the  electron  temperature  give  rough  values  of 
15  eV,  hence  the  peak  electron  and  ion  gas  dynamic 
pressure  is  approximately  4.5  x  lO^N/nr.  This 
pressure  is  higher  than  the  estimated  magnetic 
pressure(l),  1.6  x  lO^N/m^.  But,  since  it  acts 
over  a  smaller  area,  namely  the  plasma  effective 
width,  the  total  force  it  exerts  on  the  Teflon  is 
the  same  as  the  total  magnetic  force.  Recent  data (7) 
confirms  this  equal  division  when  operating  at 
the  5.5  ylb  sec  impulse  level. 


The  laser  interferometer  yields  data  up  to 
2.5  cm  from  the  Teflon  face.  The  phase  data  obtainei 
actually  gives  the  electron  density  integrated 
along  the  transverse  (z)  axis.  Density  information 
is  derived  by  assuming  that  the  transverse  variation 
is  Gaussian  with  some  estimated  effective  width. 
Figures  3  and  4  give  the  density  as  a  function 
of  time  and  of  distance  along  the  thrust  (x) 
axis.  A  two  dimensional  presentation  of  this 
data  can  be  given  by  plotting  constant  density 
contours  on  the  x-z  plane.  These  curves,  obtained 
with  the  aid  of  a  computer,  are  shown  in  Fig.  8. 

In  the  far  region,  the  time  impulsive  solutions 
apply  directly  if  the  source  duration  is  short 
compared  to  the  travel  time.  Figures  9  and  10 
show  contour  curves  of  constant  density  using  this 
solution.  In  Fig.  9  these  curves  are  connected 
smoothly  onto  the  measured  density  contours  given 
in  Fig.  8.  Since  the  far  region  distribution  function 
is  an  expanding,  drifting  Maxwellian,  the  contours 
are  expanding  circles  in  a  coordinate  system  moving 
with  the  drift  velocity  (Fig.  10). 

Several  experimental  results  agree  with  predic¬ 
tions  from  this  model.  The  theoretical  solutions 
predict  that  peak  density  should  decrease  with 
distance  by  less  than  x~l  in  the  near  region, 
by  less  than  x~2  in  the  mid  region,  and  as  x_3 
in  the  far  region.  Ion  current  is  predicted  to 
have  a  gaussian  dependence  in  the  transverse  (z) 
direction  and  an  inverse  cube  dependence  in  the 
thrust  (x)  direction.  Both  of  these  dependences 
are  verified  by  Faraday  cup  data.  Finally,  our 
model  shows  that  it  is  possible,  by  applying  convolu¬ 
tion,  to  project  downstream  densities  using  upstream 
data.  In  Fig.  7  the  density  curves  obtained  in 
this  manner  agree  in  shape  with  measured  density 
curves  and  they  differ  in  peak  magnitude  by  a 
factor  less  than  2. 

Comparing  the  Faraday  cup  data  with  the  derived 
expression  for  ion  current  in  the  far  region,  we 
estimate  a  mach  number  M  =  3,  using  the  ratio 
of  ion  drift  and  thermal  velocities.  However, 
since  the  plasma  contains  multiply  ionized  particles 
traveling  at  higher  than  average  velocities,  we 
expect  this  estimate  to  be  somewhat  high. 

The  far  region  inpul sive  solution  gives  an 
expression  for  t™,  the  time  of  peak  density  occur¬ 
rence  (see  Table  II).  For  M  =  3,  and  using  the 
value  of  tm  derived  from  the  microwave  interfer¬ 
ometer  data,  we  obtain  v0  =  27.9  km/sec,  and 
hence  Vj  *  9.3  km/sec.  These  values  agree  with 
previous  experimental  data. 

A  formula  is  derived  relating  the  total  charge 
passing  some  plane  to  the  density  at  that  plane. 
Applying  this  formula  we  obtain  estimates  of  total 
charge  on  the  order  of  that  derived  from  the  integra¬ 
tion  of  Faraday  cup  current.  Two  observations 
suggest  that  significant  recombination  takes  place 
in  the  exhaust.  First,  the  estimates  of  total 
charge  given  in  Table  III  show  a  steady  decrease 
with  distance;  and  second,  in  Fig.  7,  the  ratio 
of  peak  measured  density  to  peak  calculated  density 
decreases  with  distance. 

Questions  that  still  need  to  be  investigated 
further .  include  how  the  neutral  density  varies 
with  time,  and  what  percentage  of  multiply  ionized 
particles  are  in  the  exhaust.  These  data  are  needed 
to  improve  the  estimates  made  in  this  paper. 
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Table  II  Theoretical  Solutions  to  Time  Impulsive  Sources 


Table  III  Estimate  of  Total  Charge  Released 


Distance  from  Teflon  Face 

ij,*  (ysec) 

dg(an) 

de  (cm) 

nInt(xl016/cc) 

Q(xl0'3  Coul) 

3 . 5  mm 

1 

2.5 

.30 

6.0 

21.6 

6.0  mm 

1 

2.5 

.36 

3.1 

13.4 

11.0  mm 

1 

2.8 

.66 

1.4 

12.4 

16.0  mm 

1 

3.0 

.96 

.7 

9.7 

15.2  cm** 

5.3 

* 

As  estimated  from  Fig.  11. 

**  f  5! 

From  Faraday  Cup  Measurements K  J 
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Fig.  1(a)  Schematic  of  M-Z  interferometer  with 
source  and  photodetector. 
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lag.  1(b)  Expanded  view  of  the  recombined  beam 
showing  size  of  pinhole  in  relation 
to  period  of  intensity  variation. 
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Fig.  4  Curves  of  measured  density  as  a  function 
of  distance  along  thrust  axis  at  various 
times. 
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Fig.  5  Variation  of  maxinum  measured  density 
as  a  function  of  downstream  position. 


Fig.  7  Conparison  of  convolved  and  measured 
density  curves. 


Fig.  6  (a)  Measured  density  at  6  mm  from 

Teflon  face;  nQ(t) 

(b)  Near  region  inpulse  response; 
ni  (t) 

(c)  , (a)  Mid  region  inpulse  response; 

njft) 
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Fig.  8  Computer  plotted  constant  density 
contours  using  collected  data. 
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Fig.  9  Constant  density  contours  connecting 
data  to  theoretical  solution. 


Fig.  10  Constant  density  contours  using 
far  region  impulse  solution. 
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